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A search for long-lived charged particles reaching the muon spectrometer is performed using a data sample of 37 pb _1 
from pp collisions at y/s = 7 TeV collected by the ATLAS detector at the LHC in 2010. No excess is observed above 
the estimated background. Stable f sleptons are excluded at 95% CL up to a mass of 136 GeV, in GMSB models with 
V5 = 3, Wmessenger = 250 TeV, sign(/i) = 1 and tan/? = 5. Elcctrowcak production of sleptons is excluded up to a mass of 
110 GeV. Gluino i?-hadrons in a generic interaction model are excluded up to masses of 530 GeV to 544 GeV depending 
on the fraction of i?-hadrons produced as <?-balls. 
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1. Introduction 

Heavy long-lived particles (LLPs), with decay lengths 
longer than tens of meters, are predicted in a range of 
theories which extend the Standard Model. Supersym- 
metry (SUSY) models allow for meta-stable sleptons (I), 
squarks (q) and gauginos. Heavy LLPs produced at the 
Large Hadron Collider (LHC) could travel with velocity 
significantly lower than the speed of light. These parti- 
cles can be identified and their mass, m, determined from 
their velocity, /3, and momentum, p, using the relation 
m = p/j/3. Two different searches are presented in this 
paper, both use time of flight to measure /3, and are opti- 
mized for the somewhat different experimental signatures 
of sleptons and i?-hadrons. 

Long-lived sleptons would interact like heavy muons, 
releasing energy by ionization as they pass through the 
ATLAS detector. A search for long-lived sleptons iden- 
tified in both the inner detector (ID) and in the muon 
spectrometer (MS) is therefore performed. The results 
are interpreted in the framework of gauge-mediated SUSY 
breaking (GMSB) [l| with the light f as the LLP. If the 
mass difference between the other light sleptons and the 
light f is very small, they may also be long-lived, otherwise 
the other light sleptons decay to the f. 

Coloured LLPs (q and g) would hadronize forming R- 
hadrons, bound states composed of the LLP and light 
quarks or gluons. They may emerge as neutral states 
from the pp collision and become charged by interactions 
with the detector material, arriving as charged particles 
in the MS. A dedicated search for i?-hadrons is performed 
in which candidates are required to have MS signals while 
ID and calorimeter signals are used if available. The abil- 
ity to find i?-hadrons without requiring an ID track makes 



this analysis complementary to the previous ATLAS pa- 
per searching for i?-hadrons Q , that was based on ID and 
calorimeter signals without any requirement on the MS. 
In particular, the MS-based search presented here is more 
sensitive to models with larger <?-ball fractions. Although 
the g-ball fraction is expected to be small Q we scan the 
full range in our analysis. The results of this analysis are 
interpreted in the framework of split SUSY Q with the g 
as the LLP. 



2. Data and simulated samples 

The work presented in this paper is based on 37 pb _1 
of pp collision data collected in 2010. The events were 
selected online by muon triggers. Monte Carlo Z — > /x/x 
samples are used for resolution studies. Monte Carlo sig- 
nal samples are used to study the expected signal behav- 
ior and to set limits. The GMSB samples were generated 
with the following model parameters: the number of super- 
multiplets in the messenger sector, N5 = 3, the messenger 
mass scale, m mesS engcr = 250 TeV, the sign of the Higgsino 
mass parameter, sign(/x) = 1 and the two Higgs doublets 
vacuum expectation values ratio, tan/3 = 5. The SUSY 
particle mass scale values, A, vary from 30 to 50 TeV and 
the corresponding light f masses from 101.9 to 160.7 GeV. 
The C grav parameter was set to 5000 to ensure that the 
NLSP does not decay in the detector. The mass spectra 
of the GMSB models were generated by the SPICE pro- 
gram [B[ and the events were generated using Herwig Q. 
The i?-hadron samples were generated with g masses from 
300 to 700 GeV. As discussed in Ref. 0] several scattering 
and hadronization models can be used to describe the g R- 
hadron spectrum and interactions with the detector mate- 
rial. Three different scattering models, the first described 
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in Ref. 0, the second in Ref. [§] and the third in Ref. [§], 
and three different g-ball fractions (0.1, 0.5 and 1.0) are 
studied in this analysis. The different scattering models 
produce different fractions of candidates that arrive at the 
MS as charged particles while the <?-ball fraction affects 
the number of candidates interacting as charged particles 
in the ID. All Monte Carlo events passed the full ATLAS 
detector simulation 1(J 11 1 and were reconstructed with 
the same programs as the data. All signal Monte Carlo 
samples are normalized to the integrated luminosity of the 
data, using cross-sections calculated to next to leading or- 
der, using the PROSPINO program 
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3. The ATLAS detector 

The ATLAS detector [l3j is a multipurpose particle 
physics apparatus with a forward-backward symmetric 
cylindrical geometry and near 4n coverage in solid an- 
gle B The ID consists of a silicon pixel detector, a sili- 
con microstrip detector, and a transition radiation tracker. 
The ID is surrounded by a thin superconducting solenoid 
providing a 2 T magnetic field, and by high-granularity 
liquid- argon sampling electromagnetic calorimeters. An 
iron scintillator tile calorimeter provides hadronic cover- 
age in the central rapidity range. The end-cap and forward 
regions are instrumented with liquid-argon calorimetry for 
both electromagnetic and hadronic measurements. The 
MS surrounds the calorimeters and consists of three large 
superconducting air-core toroids each with eight coils, a 
system of precision tracking chambers, and detectors for 
triggering. 

ATLAS has a trigger system to reduce the data tak- 
ing rate from 40 MHz to ^200 Hz, designed to keep the 
events that are potentially the most interesting. The first- 
level trigger (level- 1) selection is carried out by custom 
hardware and identifies detector regions and a bunch cross- 
ing for which a trigger element was found. The high-level 
trigger is performed by dedicated software, seeded by data 
acquired from the bunch crossing and regions found at 
level- 1. The components of particular importance to this 
analysis are described in more detail below. 

3.1. The muon detectors 

The MS forms the outer part of the ATLAS detec- 
tor and is designed to detect charged particles exiting the 
barrel and end-cap calorimeters and to measure their mo- 
menta in the pseudorapidity range \rj\ < 2.7. It is also de- 
signed to trigger on these particles in the region \rj\ < 2.4. 



1 ATLAS uses a right-handed coordinate system with its origin 
at the nominal interaction point in the centre of the detector and 
the 2-axis coinciding with the axis of the beam pipe. The x-axis 
points from the interaction point to the centre of the LHC ring, and 
the y-axis points upward. Cylindrical coordinates (r, <j>) are used in 
the transverse plane, <j> being the azimuthal angle around the beam 
pipe. The pseudorapidity is defined in terms of the polar angle 9 as 
rj = - lntan(6>/2). 



The chambers in the barrel are arranged in three concen- 
tric cylindrical shells around the beam axis at radii of ap- 
proximately 5 m, 7.5 m, and 10 m. In the two end-cap 
regions, muon chambers form large wheels, perpendicular 
to the z-axis and located at distances of \z\ — 7.4 m, 10.8 
m, 14 m, and 21.5m from the interaction point. 

The precision momentum measurement is performed 
by Monitored Drift Tube (MDT) chambers, using the r\ co- 
ordinate. These chambers consist of three to eight layers of 
drift tubes and achieve an average resolution of 80 /im per 
tube. In the forward region (2 < \r\\ < 2.7), Cathode-Strip 
Chambers are used in the innermost wheel. A system of 
fast trigger chambers, consisting of Resistive Plate Cham- 
bers (RPC) in the barrel region (\rj\ < 1.05), and Thin 
Gap Chambers in the end-cap (1.05 < \rj\ < 2.4), delivers 
track information within a few tens of nanoseconds after 
the passage of the particle. The trigger chambers measure 
both coordinates of the track, rj and <j>. 

When a charged particle passes through an MDT tube 
the electrons released by ionisation drift toward the wire. 
The hit radius is obtained from the hit time, using a known 
relation between the drift distance and the drift time. A 
segment is reconstructed as a line which is tangential to the 
cylinders of constant drift distance in the different layers. 
The drift time is estimated by subtracting the muon time 
of flight, to, from the measured signal time. Slow particles 
arrive at the MDT later than muons, and if this longer 
time of flight is not taken into account, the drift distances 
are overestimated. 

The RPC chambers have an intrinsic time resolution of 
~1 ns while the digitized signal is sampled with a 3.12 ns 
granularity, allowing a measurement of the time of flight. 
When a charged particle passes through an RPC chamber 
the hit time and position are measured in the rj and <j) 
directions separately 

3.2. The tile calorimeter 

The tile calorimeter is a sampling calorimeter covering 
the barrel part of the hadronic calorimetry in ATLAS. It is 
situated in the region 2.3 < r < 4.3 m, covering \rj\ < 1.7, 
and uses iron as the passive material and plastic scintilla- 
tors as active layers. Cells are grouped radially in three 
layers. The tile calorimeter provides a timing resolution 
of 1-2 ns per cell for energy deposits typical of minimum- 
ionising particles (MIPs). The time measurement is de- 
scribed in detail in Ref. [14| . The time of flight and hence 
the velocity of a candidate can be deduced from time mea- 
surements in the tile calorimeter cells along its trajectory. 
In this analysis, only cells with a measured energy deposi- 
tion greater than 500 MeV are considered. The resolution 
of time measurements improves with increasing deposited 
energy. 
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4. Reconstruction of Long-Lived Charged Particle 
Candidates 

Penetrating LLPs leave signals similar to muons ex- 
cept for their timing, and therefore their reconstruction is 
based on muon reconstruction. However, a late-arriving 
particle may be lost in standard muon reconstruction if its 
signals are not associated in time with the collision bunch 
crossing. Late arrival of the particle also spoils segment 
fitting in the MDTs. 

The slepton search uses a dedicated muon identifica- 
tion package which starts from ID tracks and looks 
for corresponding hits in the MS, identifying candidates 
even when the segment reconstruction is imperfect, and 
refits the ID and MS hits in a combined track. Trigger de- 
tector hits arriving late with respect to the collision bunch 
crossing are also used. The next part of the LLP recon- 
struction is to estimate the particle velocity from the RPC, 
tile calorimeter and MDT [16[. The track is refitted after 
/3 has been determined, resulting in a better momentum 
resolution since it uses a set of hits which are corrected 
to take into account the late arrival of the LLP at the 
different sub-detectors. 

The i?-hadron search employs a reconstruction method 
that relies only on the MS. The reconstruction is seeded 
by a feature found by the muon trigger, without requiring 
a match with the ID. This branch of the reconstruction 
collects hits and makes segments starting from the posi- 
tion and momentum of the trigger candidate in the mid- 
dle station of the MS and extrapolates the track to the 
other stations. Once all segments are reconstructed, /3 is 
estimated. A candidate which is not found by the muon 
trigger, i.e. if it arrives late at the trigger chambers and its 
hits are not associated with the collision bunch crossing, is 
not reconstructed in the MS-standalone method, leading 
to a loss of efficiency at low (3. 

4-1. j3 estimation 

The value of j3 for each candidate is estimated by mini- 
mizing the total \ 2 between the available timing measure- 
ments from the sub-detectors, and the timing expected 
from the hypothesized (3 value. Contributions to the \ 2 
are calculated as follows. 

MDT segments: An MDT segment is reconstructed 
as a line tangent to the circles of constant drift distance 
in the different layers, after the radii are estimated from 
the drift time using a known relation R (Adrift), where tdrift 
is estimated as ^measured — to- Slow particles have a longer 
time of flight, and a better segment fit is obtained with the 
correct to. For each test (3, new MDT segments are built 
from a set of hits in a road around the extrapolated track, 
using the to corresponding to the arrival time of a particle 
traveling with velocity /3, R{t mclisurc a - t (/3)). The \ 2 
representing the difference between the inferred position 
of the particle in the set of tubes and the segment position 
in the tube is minimized. For a low f3 particle, this method 
recovers hits on segments that could be lost when fitting 



the segment assuming /3 = 1. This estimate of the MDT 
f3 is used in the i?-hadron search. 

MDT hits on track: For the slepton search, the time 
estimate from the MDT segment method can be improved 
by performing a full track fit to the ID and MS hits. The 
estimated particle trajectory through each tube is signif- 
icantly more accurate after the full track fit than in the 
segment finding stage. The time of flight of the parti- 
cle to each tube is obtained using the difference between 
the time of flight corresponding to the refitted track po- 
sition in the tube, £r and the time actually measured, 
to = ^measured — £r- The \ 2 between the measured time 
of flight and the time of flight corresponding to the arrival 
time of a particle traveling with the test {3 is minimized. 

RPC and tile calorimeter: The position and time 
are independent measurements for each hit (or cell). The 
X 2 minimization is performed using the measured times of 
hits on the candidate track. 

The time of flight measurement quality is sensitive to 
the time resolution of the detector. In a perfectly cali- 
brated detector, any energetic muon coming from a col- 
lision in the interaction point will pass the detector at 
to = 0. The to distributions in the different sub-detectors 
are measured and their means used to correct the cali- 
bration. The observed width of these distributions after 
correction is used as the error on the time measurement in 
the (3 fit and to smear times in the simulated samples. 

The (3 distributions of candidates obtained in the com- 
bined minimization are shown in Figure [TJ after the (3- 
quality selection described in section [5] For the slepton 
search the mean [3 value is 0.997 and the resolution is 
o-p = 0.048. For the _R-hadron search the mean value is 
(3 = 1. 001 and the resolution is ap — 0.051. 

4-2. Signal resolution expected in data 

Since (3 is estimated from the measured time of flight, 
for a given resolution on the time measurement, a slower 
particle has a better f3 resolution. To simulate correctly 
the time resolution corresponding to the current state of 
calibration, hit times in simulated samples are smeared to 
reproduce the resolution measured in the data, prior to 
the j3 estimation. Figure [T] shows the j3 distribution for 
selected Z — > fi/i decays in data and in Monte Carlo with 
smeared hit times. It can be seen that the smearing mech- 
anism reproduces the measured muon f3 distribution. The 
same time-smearing mechanism is applied to the signal 
Monte Carlo samples. 

5. Candidate selection 

5.1. Trigger selection 

This analysis is based on events collected by two types 
of muon trigger chains. The trigger for the slepton search 
requires MS tracks to be matched with ID tracks in the 
high-level trigger. The estimated pt is obtained from the 
combination of both systems, and is required to satisfy 
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Figure 1: Distribution of j3 for all candidates in data (points with 
error bars), muons from the decay Z — ¥ /ifi in data (full lines) and 
smeared Monte Carlo (dashed lines), in the estimation used in the 
slepton search (upper) and in the estimation used in the .R-hadron 
search (lower). 



Pt > 13 GeV. The trigger for the i?-hadron search requires 
an MS-standalone muon trigger with pt > 40 GeV. The 
standalone triggers have less accurate pt estimates than 
the combined ones. 

The events are selected online by requiring at least 
one level-1 muon trigger. As level-1 muon triggers are 
accepted and passed to the high-level trigger only if as- 
signed to the collision bunch crossing, late triggers due to 
late arrival of the particles are lost. However, late trig- 
gers due to late arrival of slow particles from the collision 
bunch crossing are also lost. The level-1 trigger efficiency 
for particles arriving late at the MS is difficult to assess 
from data, where the overwhelming majority of candidates 
are muons. Therefore the trigger tracking efficiency is es- 
timated from Z — > fifi data, while the effects of timing 
are estimated from simulated i?-hadron and GMSB events 
passing the level-1 trigger simulation, which includes the 
timing requirements of the trigger electronics. The es- 
timated trigger efficiencies for GMSB slepton candidates 
are between 80% and 81%. The i?-hadron search is much 
more adversely affected by the loss of trigger efficiency for 
late candidates, since the reconstruction is seeded by the 
trigger, and so even if an event is triggered by another ob- 
ject, the candidate will be lost. For i?-hadrons that could 
be reconstructed because they are charged in the MS, the 
trigger efficiencies range between 55% for nig = 300 GeV 
to 38% for nig = 700 GeV. The estimated trigger efficiency 
with respect to all i?-hadrons produced in the scattering 
model of Ref. [7J varies from 25% for m g = 300 GeV to 
17% for nig = 700 GeV. The effect of the trigger efficien- 
cies can be seen in Tables [1] and [5] for the signal and data 
in the slepton and i?-hadron searches respectively. 

5.2. Offline selection 

Collision events are selected by requiring a good pri- 
mary vertex with more than two ID tracks, and with |zj tx | < 
150 mm (where Zq' x is the z coordinate of the recon- 
structed primary vertex) . 

Cosmic ray background is rejected by removing tracks 
that do not pass close to the primary vertex in z. For 
candidates with an associated ID track, candidates with 
| ^o rk — z o tx | > 10 mm are removed, where Zy k is the z co- 
ordinate at the distance of closest approach of the track to 
the origin. If no ID track is associated with the candidate, 
then it is still rejected if |zQ lk — Zg tx | > 150 mm. Pairs of 
candidates with approximately opposite n and <j> are also 
removed. 

The analysis searching for sleptons requires two can- 
didates in each event, because two sleptons are produced, 
and both have a high probability to be observed in the 
MS. However, only one of them is required to pass the LLP 
selection. This requirement reduces background from W 
production and QCD, but Z — > fifi decays remain. Any 
candidate that, when combined with a second muon, gives 
an invariant mass within 10 GeV of the Z mass is rejected. 
In the i?-hadron search, no requirement of two candidates 
per event is made, because i?-hadrons may be neutral in 
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the MS, or be lost by triggering in the next bunch cross- 
ing. Nevertheless, pairs consistent with the Z mass are 
still rejected in the i?-hadron search. The above require- 
ments are grouped in Tables Q] and [2] under the label "event 
selection" . 

The slepton search requires candidates to have pt > 
40 GeV, well above the efficiency plateau for the trigger 
threshold of 13 GeV. A pt requirement of 60 GeV is ap- 
plied for all candidates in the iZ-hadron search, so as to 
be in the MS-standalone trigger-efficiency plateau. Can- 
didates with pt > 1 TeV are rejected. This removes a 
few candidates with badly reconstructed momenta in both 
searches. Each candidate is required to have \t)\ < 2.5. 
These requirements are grouped in Tables [T] and [2] under 
the label "candidate quality" . 

The estimated ft is required to be consistent for mea- 
surements in the same sub-detector, based on the RMS 
of ft calculated from each hit separately. The estimated 
ft is also required to be consistent between sub-detectors. 
A ft measurement in at least two subdetectors is required 
for \r)\ < 1.7. These requirements are grouped in Tables [1] 
and [2] under the label 11 ft quality" . Finally, in order to re- 
ject most muons, the combined ft measurement is required 
to be in the range ft < 0.95. 

6. Background estimation 

The background is mainly composed of high muons 
with mis-measured ft. The estimation of the background 
mass distribution is made directly from the data and re- 
lies on two premises: that the signal to background ratio 
before applying requirements on ft is small and that the 
probability density function (pdf) for the ft resolution for 
muons is independent of the source of the muon and its 
momentum. 

For each muon candidate passing the ft quality require- 
ment, a random ft is drawn from the muon ft pdf. If this 
ft is inside the signal range, j3 < 0.95, a mass is calculated 
using the reconstructed momentum of the muon and the 
random ft. The statistical error of the background estima- 
tion is reduced by repeating the procedure many times for 
each muon and dividing the resulting distribution by the 
number of repetitions. The mass histogram obtained this 
way represents the background estimation. 

The ft distribution is different in different detector re- 
gions for three main reasons: different r\ regions are cov- 
ered by different technologies (\r]\ < 1.05 for RPC, \rj\ < 
1.7 for tile calorimeter, \r]\ < 2.5 for MDT); the time of 
flight method is more precise when the distance between 
the interaction point and the detector element is larger; 
the measurement in some regions of the detector is less 
precise due to fewer detector layers and magnetic field in- 
homogeneities. The background estimation is performed 
in n regions so that the ft resolution within each region 
is approximately the same. The muon ft pdf in each rj 
region is given by the histogram of the measured ft of all 
muons in the region. The regions also differ in the muon 



momentum distribution, since for any cut, p is larger 
as r) increases. Therefore the combination of p with ran- 
dom ft is done separately in each region and the resulting 
distributions are added together. 

The efficiency of the requirements described in sec- 
tion 15.21 to reject cosmic rays is estimated from data col- 
lected with a cosmic muon trigger in the empty bunches 
and periods without collisions, dropping the requirement 
of a good primary vertex. The number of remaining cosmic 
ray muons are estimated from the number of candidates 
rejected by these requirements in the collision sample and 
the rejection efficiency. This results in 1.3±0.2 cosmic rays 
in the i?-hadron search sample. The estimated cosmic ray 
contamination in the slepton search sample is 0.7 ± 0.2 
candidates. 

7. Systematic uncertainties 

Several possible sources of systematic uncertainty have 
been evaluated. 

7.1. Signal yields 

The total experimental systematic uncertainty in the 
signal yields is 6% on average. The sources and their in- 
dividual contributions are described below. 

An uncertainty of 3.4% is assigned to the measurement 
of the integrated luminosity [l7j . 

The systematic uncertainty associated with the trigger 
selection is estimated in Ref. [U^ to be 0.73% (0.35%) and 
0.74% (0.42%) in the barrel (endcap) for the two trigger 
chains used in the slepton search. An uncertainty of 5% is 
estimated for the i?-hadron search using similar methods. 
The trigger simulation models in detail the timing require- 
ments of the trigger electronics and the efficiency loss due 
to particles arriving late at the MS. Differences between 
data and Monte Carlo trigger efficiency due to the differ- 
ence in time resolutions between data and simulation were 
tested and are negligible. 

The signal ft resolution expected in the data is esti- 
mated by smearing the hit times according to the spread 
observed in the time calibration. The systematic uncer- 
tainty due to the smearing process is estimated by scal- 
ing the smearing factor up and down, so as to bracket 
the distribution obtained in data. A 6% (2%) system- 
atic uncertainty is associated with the smearing process 
for the GMSB models in the barrel (endcap). For the g 
i?-hadrons, the effect of the smearing is negligible. 

The systematic uncertainties due to track reconstruc- 
tion efficiency and momentum resolution differences be- 
tween ATLAS data and simulation are estimated to be 
0.5% for GMSB events and between 0.8% and 1.3% for 
i?-hadrons in the different hadronization and interaction 
models. 
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Table 1: Candidates in data compared to the simulated GMSB signal passing the selection stages in the slepton search. The Monte Carlo 
signal prediction is normalized to the data luminosity using the next to leading order cross-section. 
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Table 2: Candidates in data compared to the simulated ij-hadron signal passing the selection stages in the ij-hadron search. The Monte 
Carlo signal prediction for the sample with the scattering model of Ref. Q and a g-ball fraction of 0.1 is normalized to the data luminosity. 



7.2. Background estimate 

A total of 15% (20%) uncertainty on the background is 
estimated for the slepton (i?-hadron) search resulting from 
individual contributions discussed below. 

A systematic variation of the /3 distribution within each 
of the detector regions used in the background estimation 
may lead to a systematic error on the background esti- 
mation. To quantify the variability of the f3 distribution 
within a region and its effect on the background estimate, 
each region is sub-divided into smaller regions and the 
variation of their /3 distribution is used as a variability 
estimate. This leads to the dominant uncertainty in the 
background estimate. 

The possibility that a /3-distribution dependence on the 
candidate momentum or source would result in a system- 
atic uncertainty on the background estimate was tested. 
The candidates in each r\ region were divided by their mo- 
mentum into two bins with similar counts. The indepen- 
dence of the j3 pdf from the source of the muons is con- 
firmed using Z —¥ fifi samples. Estimating the background 
using the pdf from the low or high momentum bins or from 
muons from Z — > fifi results in negligible systematic un- 
certainties. 

Finally, the background estimation is based on a lim- 
ited statistics sample, that of all candidates that passed 
the candidate quality requirements, before the cut on j3. 
The tail of the background mass distribution has a sig- 
nificant contribution from a few high momentum events, 
and a statistical error arises from this. In order to cal- 
culate the sensitivity of the limits to the statistics of the 
momentum distribution, the candidate sample was divided 
randomly into two samples and the background estimate 



derived from each sample separately. The resulting uncer- 
tainty in the slepton search ranges from 1.04 candidates for 
m > 90 GeV to 0.14 candidates for m > 140 GeV, while 
the errors from sample statistics in the i?-hadron search 
are negligible. 

7.3. Theoretical cross-sections 



The PROSPINO program [12J is used to calculate the 
signal cross-sections at next to leading order and two sources 
of theoretical systematic uncertainties were considered: the 
renormalization and factorization scales are changed up- 
ward and downward by a factor of two. This results in a 
systematic error of 7% for GMSB cross-sections and 15% 
for g cross-sections [2|. The parton density functions of 
CTEQ6.6 1!) are used, and the uncertainty due to varia- 
tions in the parton distribution functions is estimated to 
be 5%. 

8. Results 

Figure [5] shows the candidate mass distribution for 
data and the estimated background with its systematic 
uncertainty. Good agreement is observed. The CL S ap- 
proach 20( for counting experiments is used to derive the 
limits for the production cross-section of GMSB slepton 
and g i?-hadron events. The limits are obtained by com- 
paring the expected number of events with a candidate 
above a given mass cut with the actual number of events 
with a candidate above the same mass cut observed in the 
data. For each model, the mass cut is chosen to give the 
best expected limit. The mass cuts for the different mod- 
els are summarized in Tables [3] and @] together with the 
expected signal and background in each case. 
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Table 3: Mass cut and expected number of events as a function of the 
f mass in the slepton search. The systematic uncertainties on the 
signal yield and background estimate are 6% and 15% respectively. 
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Table 4: Mass cut and expected number of events as a function of the 
g mass in the ij-hadron search. The systematic uncertainties on the 
signal yield and background estimate are 6% and 20% respectively. 



The expected number of signal candidates for an inte- 
grated luminosity of 37 pb -1 is added to the background 
estimation and compared to the data in Figure [2] for the 
GMSB and i?-hadron models. The production cross-section 
for f events and cross-section limit is shown in Figure [3J 
as a function of the f mass for the slepton search (upper). 
Stable f are excluded at 95% CL up to a mass of 136 GeV, 
in GMSB models with N$ — 3, m mesS enger = 250 TeV, 
sign(/i) = 1 and tan/3 = 5. Figure [3] (lower) shows the 
limit obtained for sleptons produced only by electroweak 
processes, which have a smaller dependence on the model 
parameters other than the slepton mass. Sleptons pro- 
duced in electroweak processes are excluded up to a mass 
of 110 GeV. Previous limits on stable sleptons are all be- 
These limits are only applicable to mod- 



low 100 GeV 21 



els where the f or sleptons are the next to lightest SUSY 
particle, and their lifetime is sufficiently long to traverse 
the ATLAS experiment. In this case, the limits obtained 
for the above models are expected to have limited depen- 
dence on tan/3 and N$. 

Figure @] shows the limits for g i?-hadrons in the scat- 
tering model of Ref. Q ■ Such i?-hadrons are excluded at 
95% CL up to a mass of 544 GeV for a g-ball fraction of 
0.1. Models with <?-ball fractions of 0.5 and 1.0 are ex- 
cluded up to masses of 537 GeV and 530 GeV respectively. 
Previous less stringent limits were set by the Tevatron [22j | 
and by the CMS Collaboration [23] which used the MS to 
select candidates, but not to measure /3. For a <?-ball frac- 
tion of 0.1, the ATLAS collaboration in Ref. Q sets limits 
that are higher then the limits presented here, not using 
the MS. 



-♦- data 

— background estimate: 
GMSB 

101.9 GeV 
11 6.3 GeV -= 
131.0 GeV 




150 200 250 300 
m [GeV] 



data 

background estimate; 

R-Hadrons „„ „ , , 
*00 GeV 
rrfi = 400 GeV 
m~° = 500 GeV 
" 600 GeV 




100 200 300 400 500 600 700 
m [GeV] 



Figure 2: Candidate estimated mass distribution for data, expected 
background including systematic uncertainty, with simulated signals 
added, in the slepton (upper) and ij-hadron (lower) searches. 
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Figure 3: The expected production cross-section for GMSB events 
with N$ = 3, m me ssenger = 250 TeV, sign(/^) = 1 and tan/3 = 5, 
and the cross-section upper limit at 95% CL for the slepton search 
as a function of the f mass (upper) and for sleptons produced in 
electroweak processes only (lower). 
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Figure 4: The expected production cross-section for ij-hadron events 
and the cross-section limit at 95% CL as a function of the g mass 
for the i?-hadron search in the scattering model of Ref. Q and for 
different g-ball fractions. The expected limit and its 1 a band are 
shown for 0.1 g-ball fraction. 



9. Summary and conclusion 

A search for long-lived charged particles reaching the 
muon spectrometer was performed with 37 pb -1 of data 
collected with the ATLAS detector. No excess is observed 
above the estimated background and 95% CL limits on f 
and iZ-hadron production are set. Stable f's are excluded 
up to a mass of 136 GeV, in GMSB models with N$ = 3, 
"^messenger = 250 TeV, sign(/x) = 1 and tan/? = 5. Slep- 
tons produced in electroweak processes are excluded up to 
a mass of 110 GeV. Gluino i?-hadrons in the scattering 
model of Ref. Q are excluded up to masses of 530 GeV to 
544 GeV depending on the fraction of i?-hadrons produced 
as g-balls. 
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